An Infrared study of the Josephson vortex state in high-Tc cuprates 
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We report the results of the c-axis infrared spectroscopy of La2-^a:Sra:Cu04 in high magnetic field 
oriented parallel to the Cu02 planes. A significant suppression of the superfluid density with mag- 
netic field ps{II) is observed for both underdoped (x=0.125) and overdoped (x=0.17) samples. We 
show that the existing theoretical models of the Josephson vortex state fail to consistently describe 
the observed effects and discuss possible reasons for the discrepancies. 
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Although the mechanism of high-Tc superconductivity 
in cuprates still remains unresolved, significant progress 
has been made in understanding of the mixed state of 
these systems jl). Many subtle predictions of the theo- 
ries of the vortex state have been experimentally verified 
^]-p7|. In particular, electrodynamics in relatively small 
magnetic fields (typically H<7 T) appears to be in ac- 
cord with the models proposed for both pancake vortices 
formed with H ± Cu02-planes, and for the Josephson 
vortices developing in the H || Cu02-plane configuration 
In this work we explored field-induced suppres- 
sion of the superfluid density in La2_a;Sra;Cu04 (LSCO) 
single crystals due to Josephson vortices. By extend- 
ing previous studies to stronger fields (up to 17 T) we 
find marked departures of the experimental data from 
conventional Josephson vortex theory. In particular the 
superfluid density in undcrdoped LSCO appears to be 
much more fragile than the models prescribe. The impli- 
cations of these flndings for our current understanding of 
cuprates will be discussed. 

We used infrared (IR) spectroscopy to examine the in- 
terlayer (E || c) response of high quality LSCO crystals 
[|8|: an underdoped with x=0.125 (Tc ~ 32 K) and a 
weakly overdoped with x=0.17 (Tc ~ 36 K). The near- 
normal-incidence reflectance -R(w) was collected over a 
broad temperature (6-300 K) and frequency (10-48,000 
cm~^) range. In addition field-induced changes of the 
refiectivity R{H,6K)/ R{0T,6K) were measured under 
zero field cooling with H || Cu02, when vortices pen- 
etrate in-between the Cu02 planes The uncer- 
tainty of the absolute measurements in our apparatus 
does not exceed 0.5 %, whereas the relative errors of 
the field-induced changes are less than 1 %. The optical 
conductivity cri(a;) -I- i(T2(w) and the dielectric function 
ei(a;) -|-ie2(w) were calculated from R{uj) using Kramers- 
Kronig analysis. 

In the superconducting state the real part of the 
optical conductivity has two components: crf'^(ijj) ~ 
^S{0) + cr['^^((jj), the first one due to superconducting 
condensate and the second due to unpaired carriers be- 
low Tj, (a so called regular contribution) . The superfluid 
density ps is quantified with the plasma frequency ajg-. 
Ps = ujI = (? iy?f. = Ane^ns/m* related to the density of 
superconducting carriers n^ and their mass m* ; Ac is the 



penetration depth. A delta function in crj' [lo) gives rise 
to a term — (ws/w)^ in ei(w). A common procedure of 
extracting lOs from ei(ci;) involves fitting of the low en- 
ergy part of the spectrum with l/uP . The problem with 
this procedure is that it does not discriminate screening 
due to superconducting condensate from regular contri- 
bution. In order to fix this problem we use the following 
correction procedure: 



ei(w) - e™^(w) 



w2 
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where {lo) is the regular contribution to dielectric 
function. In order to calculate e"^^ {oj) a KK-like trans- 
formation is employed: 
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where e~2^{Lj) is the regular contribution to the imaginary 
part of dielectric constant, i.e. after a (5(0)-function has 
been subtracted. The procedure described by Eq. ^ and |^ 
accounts for contribution of unpaired carriers at T < Tc, 
but also phonons, interband transitions, magnons, and 
all other finite energy excitations. 

Figure ^ displays the raw reflectance data. In zero field 
the 6 K reflectance of both LSCO crystals is character- 
ized by a sharp plasma edge at the frequency lOs I \J^oa 
(too is the high frequency dielectric constant). This form 
of reflectance, which resembles a plasma edge in normal 
metals, is due to the zero crossing in the real part of the 
dielectric function, txiui) = 0, and is commonly referred 
to as the Josephson Plasma Resonance (JPR). As tem- 
perature increases, the edge in Rito) is smeared out and 
the minimum shifts to lower frequency, indicating sup- 
pression of the superfluid density. In the normal state the 
underdoped sample shows a very weak upturn of R{ljS) 
as w 0, indicating a small metallic contribution at 
T > Tc. The upturn is stronger in the x=0.17 crystal 
suggesting that the far IR conductivity increases with 
doping. Similar enhancement of the "metallic" trends 
in the interlayer response is commonly found in other 
cuprate families. 

Application of a magnetic field has a strong impact on 
the JPR feature in both crystals (bottom panels): the 
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plasma edges are smeared and the minima in R{u>, H) 
also shift to lower w. The plasma edge shift in 17 T field 
is as strong as 35 % in the x=0.125 sample and 10 % 
in the x=0.17 crystal. This result is surprising given the 
fact that the strongest field used in our experiments is 
still orders of magnitude smaller than the upper critical 
field Hc2 for the H ||Cu02 orientation. 

It is instructive to present the evolution of the super- 
conducting condensate with H and T with the spectra 
of the loss function /?7i(l/e(a;)) = e2(w)/(ei(w) + eiC"^)) 
(Fig. H). At T <C Tc the loss function is peaked at the 
frequency close to the JPR, whereas the width of the 
/m(l/e(w)) mode is proportional to the magnitude of 
e2(w,T < Tc). As temperature increases the peak shifts 
to lower energies, indicating suppression of the superfluid 
density. The insets in the top panels show the tempera- 
ture dependence of the superfluid density closely resem- 
bling the form expected for a d-wave superconductor. At 
the lowest measured temperature we found the following 
values of the superfluid density: u;s=160 cm^^ (Ac = 9.7 
^m) for 12.5 % sample and a;s=360 cm~^ (Ac = 4.3 ^m) 
for 17 % material. 

The behavior of the loss function in high magnetic field 
(Fig. H bottom panels) is generally similar to the H=0 
data taken at finite temperature: the peak softens and 
its width is enhanced as the field increases. Fig.3A quan- 
tifies the demise of lOs in magnetic field. At 17 T the 
superfluid density is reduced by 38 % in the underdoped 
and 12 % in the overdoped sample; these values are in full 
agreement with the strength of the effect inferred directly 
from raw data in Fig. ^ We also note qualitatively differ- 
ent form of the uJs{H) dependencies in the two samples. 
Two principal mechanisms of the superfluid density sup- 
pression in high magnetic field are: 1) direct pair-braking 
of the condensate and 2) dissipation associated with vor- 
tex dynamics. The former process is usually discarded in 
view of giant -ffc2 values in the H || Cu02 configuration. 
Within the latter picture the oscillating electric field with 
the E-vector along the c-axis leads to a transverse mo- 
tion u of Josephson vortices ^ located between the Cu02 
planes in the direction perpendicular to the field (see the 
bottom inset of Fig. ^). Electrodynamics of Josephson 
vortices has been thoroughly discussed in several papers. 
[||-|l6| Below we show that the behavior of JPR in LSCO 
crystals cannot be fully understood within the proposed 
models, especially in the underdoped sample. 

Bulaevskii et al have worked out the following predic- 
tion for the field dependence of uJs '■ [p6| 



1--— In^ 
8 Ho ^ H 
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The effect of magnetic field is directly related to the 
strength of the coupling between the Cu02 layers 
that is parameterized through the characteristic field 
Ho=$o/7s^- In this expression $0 is the fiux quantum, 



7 = Xc/Xab (Aa6 is the in-plane penetration depth) is the 
anisotropy factor and s ~ 1.32 nm is the interlayer dis- 
tance. We estimated Ho=36.7 T for the x=0.125 with 
Xab = 0.2 /xm (Ref. |l|) and Ho= 55 T for the x=0.17 
sample with Aob = 0.3 ^m (Ref. |l^]). In Fig. 3B we 
plot uJs(H)/ijJs{0) as a function of H/Hg along with the 
theoretical dependence (full line). The overall trends of 
oJs{H) are clearly different for the two crystals. In partic- 
ular, the x=0.125 material reveals a dramatic reduction 
of the condensate strength compared to the model pre- 
diction for H/Hq > 0.3. For x=0.17 sample such high 
magnetic fields could not be achieved with the present 
experimental setup. 

We also attempted to describe uJs{H) within a 
phenomenological scenario of vortex dynamics |]l^ , |l5| . 
Tachiki, Koyama and Takahashi (TKT) ||l^ obtained an 
explicit result for the complex dielectric function of a lay- 
ered superconductor in the mixed state: 



<t>a 
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where M is the vortex inertial mass, 77 is the viscous 
force coefficient and tip IS the vortex pinning constant. 
In Eq. ^, is the real part of the dielectric function 
above the plasmon and a;„ and jsr are the regular com- 
ponent (due to unpaired carriers below Tc) plasma fre- 
quency and scattering rate. Eq. ^ can be regarded as a 
H^O generalization of the well known " two fluid" model 
of superconductivity, commonly used in the microwave 
and IR frequency ranges |^ . In order to model the data 
in Fig. |3| we extracted the ujs{H) behavior from Eq. |4|by 



exploring Re[e(ci;)] in the limit of — > 0: lu^ 



'xei(cj). 



Some of the parameters of the TKT equation needed to 
obtain this result were readily available from the fits of 
R{lu,6K,H = 0): €00=27, w„=200cm-i (for 12.5 % sam- 
ple) and 1,100 cm^^ (for 17 % sample), jsr— 5,000 cm^^. 
Similar to all previous spectroscopic works we set 

the vortex mass to zero in Eq. ^. The viscous drag con- 
stant 77 can be calculated within Coffey-Clem approach 
yielding ri—7 Pa cm and 28 Pa cm, for the under- 
doped and overdoped samples respectively. There- 
fore we are left with the pinning force constant f^p clS 
the only fltting parameter in Eq. |[ As Fig. 3A shows 
the theory gives a very good flt for the overdoped sam- 
ple, with Kp= 6,000-11,000 Pa, a value comparable with 
that of nearly optimally doped YBa2Cu30j, (in the 
same field configuration). For the underdoped sample on 
the other hand, we could not obtain a good fit for any 
value of Kp . In order to reproduce the overall depression 
of the superfluid density in 17 T field we have to adopt 
Kp=150-200 Pa. We believe that such a vast difference in 
the magnitude of Kp between the two samples is another 
signal of inability of the TKT scenario to account for 
the experimental situation at least in underdoped LSCO 
crystals. 
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One obvious distinction between the underdoped and 
overdoped samples is the different nature of the inter- 
layer transport. Earlier experiments |Q established that 
the x=0.16 doping in the LSCO system separates the re- 
gion of " insulating" (forx<0.16) and "metallic" (x>0.16) 
ground states. This result may have important impli- 
cations for the nature of the H || Cu02-vortices. The 
x=0.125 crystal is likely to fall in the regime where a 
description within the formalism of Joscphson vortices 
is applicable owing to the insulating character of the 
"medium" separating the Cu02 planes. However, as dop- 
ing progresses to the over-doped side eventually the H | 
Cu02-vortices will evolve into Abrikosov-type vortices 
having a normal core. Therefore our x=0.17 crystal may 
be closer to the regime where additional factors involving 
complex character of the vortex cores in d-wave supercon- 
ductors have to be taken into consideration. |p^ , p5| , p6| 
Surprisingly, the Josephson analysis (Eqs. js] -^) obvi- 
ously ignoring the latter issues, is less problematic for the 
x=0.17 crystal but fails for the x=0.125 sample which, 
based on its "insulating" resistivity at T 0, is like- 
lier to comply with Eqs. ^ - In particular, the TKT 
analysis yields an unexpectedly small pinning constant 
which is hard to reconcile with prominence of the intrin- 
sic inhomogeneities in underdoped cuprates, since inho- 
mogeneities would normally promote vortex pininnig. 

When searching for reasons for the inability of Eqs. ^ 
- ^ to consistently describe the anomalous sensitivity of 
the superfluid response in underdoped LSCO to mag- 
netic field it is prudent to revisit the assumptions of these 
models. Indeed, Eqs. § - ^ are valid for Josephson cou- 
pling between uniform s-wave superconductors; validity 
of both assumptions for high-Tc materials is question- 
able. One implication of the c?-wave order parameter 
(firmly established for cuprates) is that the Zeeman en- 
ergy associated with the H || Cu02 field can no longer 
be disregarded for states close to the node. The im- 
pact of the Zeeman field is especially important for the 
x=0.125 phase since at this particular composition the 
magnitude of the in-plane superconducting gap is anoma- 
lously low. The field-dependence of the c-axis super- 
fluid density within this scenario has the following form: 
w2(i/)/cj2(o) = ^1 - (H/A)2 [H and is displayed in 
Fig.SC for different magnitudes of the energy gap 2A. 
Interestingly, the functional form of lOs{H) / uOs{H — 0) is 
reproduced by this calculation. The Zeeman reduction of 
the condensate strength occurs in concert with other dis- 
sipation mechanisms and accounts for at least 10-15% de- 
pression of the a;s(17T) value. An additional factor per- 
tinent to the anomalous field response of the underdoped 
crystals may be connected to spatial non-uniformities of 
superconductivity within the Cu02 planes revealed by a 
variety of experimental methods. pO| . It seems plausible 
that the magnetic field will infiuence coupling between 
these dissimilar regions in the Cu02 plane thus leading 
to field-dependence of the in-plane superconducting pa- 



rameters such as Xab in Eq. ^. 

Models of the superfluid density suppression discussed 
above involve only a reduction of carrier density n^ . How- 
ever, the c-axis response of cuprates also reveals changes 
of the effective mass m* or of the kinetic energy at 
T <Tc. Recently loffe and Millis proposed that the 
variation the effective mass below Tc is connected with 
phase coherence between the Cu02 planes. |Q Phase co- 
herence is suppressed in magnetic field which may lead 
to a more rapid degradation of the superfluid density via 
the increase of m* in addition to usual reduction of n^. 
This latter effect is expected to be particularly strong in 
underdoped samples which show the strongest changes of 
TO* at T < Tc. 

In conclusion, magneto-optical results for underdoped 
LSCO revealed remarkable depression of the superfluid 
density in the vortex state. We have identified sev- 
eral factors which may account for much more com- 
plex behavior of underdoped LSCO beyond conventional 
Josephson vortex models. Further theoretical analysis is 
needed to distinguish between the roles played by spatial 
non-uniformities of superconducting state, as well as by 
changes of kinetic energy in the observed behavior. In 
this fashion quantitative understanding of vortex state 
data presented in this work will be instrumental in nar- 
rowing down the range of plausible theoretical descrip- 
tions of the miderdoped state cuprates. 
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FIG. 1. Infrared data for LSCO crystals with x=0.125 (left 
panels) and x=0.17 (right panels). The top panels show c-axis 
reflectance in zero field; the bottom panels: R{uo) in high 
magnetic field. 

FIG. 2. The loss function Im(l/e(cj)) reveals coupling to 
the longitudinal JPR mode. The top two panels show the 
temperature dependence of the loss function and the bottom 
ones the loss function in high magnetic field. The top insets 
show the temperature dependence of the superfluid density 
and the bottom one sketch the geometry of a Josephson vortex 
experiment. 

FIG. 3. Change of the superfluid density with magnetic 
field uos[H) /ujsiO) alone with the theoretical results obtained 
for the TKT model (panel A), Bulaevskii et al model (panel 
B) and scenario taking into account the nodal Zeeman effect 
(panel C). While conventional models of the Josephson vortex 
state are inconsistent with the high field data for the under- 
doped x=0.125 crystal (A and B), plausible description can be 
achieved within the picture discussed by Won et al. [^ij using 
the magnitude of the gap in from the in-plane measurements 
for the same crystal [BSl. 
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